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By Claxence A. Syvertson and David H. Dennis 


SUMMARY 


A second-order shock-expansion method applicable to bodies of revolu- 
tion near zero lift is developed. Expressions defining the pressures on 
noninclined bodies are derived by the use of characteristics theory in 
combination with properties of the flow predicted by the generalized 
shock-expansion method. This result is extended to inclined bodies to 
obtain expressions for the normal-force and pitching-moment derivatives 
at zero angle of attack. The method is intended for application \mder 
conditions between the ranges of applicability of the second-order poten- 
tial theory and the generalized shock-expansion method - namely, when the 
ratio of free-stream Mach number to nose fineness ratio is in the neighbor- 
hood of 1. 

For noninclined bodies, the pressure distributions predicted by the 
second-order shock-expansion method are compared with existing experimental 
resijlts and with predictions of other theories. For inclined bodies, the 
normal-force derivatives and locations of the center of pressiu'e at zero 
angle of attack predicted by the method are compared with experimental 
results for Mach, numbers from 3.00 to 6.28. Fineness ratio 7, 5, and 3 
cones and tangent ogives were tested alone and with cylindrical afterbodies 
up to 10 diameters long. In general, the predictions of the present method 
are found to be in good agreement with the experimental restilts. For non- 
inclined bodies, pressiore distributions predicted with the method are in 
good agreement with existing experimental resvilts and with distributions 
obtained with the method of characteristics. For inclined bodies, the 
normal-force derivatives per radian (for normal-force coefficients refer- 
enced to body base area) are predicted within ±0.2 and the locations of 
the center of pressure are predicted within ±0.2 body diameters. On the 
basis of these results, the second-order shock-expansion method appears 
applicable for values of the ratio of free-stream Mach number to nose 
fineness ratio from 0.4 to 2. 


INTRODUCTION 


The flow about bodies traveling at high supersonic speeds was inves- 
tigated by Eggers (ref. 1). He found that under specified conditions such 
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flows could be considered as locally two-dimensional and that they could 
be treated by a generalized shock- expans ion method. The application of 
this method to nonlifting bodies of revolution had been given previously 
(ref. 2), and subsequently the method was applied to lifting bodies in 
references 3 It was found that the generalized shock-expansion 

method accurately predicted the flow about pointed bodies of revolution 
when the hypersonic similarity parameter (ratio of Mach number to body 
fineness ratio) was greater than about 1. This method is, therefore, 
particiilarly useful in the treatment of flows about bodies traveling at 
relatively lar'ge Mach nvanbers . At lower speeds , the second-order potential, 
theory of Van Dyke (ref. 5) has been widely used. (See, also, his hybrid 
theory for slightly inclined bodies, ref. 6.) The application of this 
theory to bodies traveling at large Mach nimibers is often limited, however, 
by the restriction that the maximum slope of the body must be somewhat less 
than the slope of a free-stream Mach line. 

The ranges of applicability of the generalized shock-expansion method 
and the second-order potential theory do not always overlap, and there 
remain, therefore, flows at certain combinations of Mach number and body 
shape which cannot be treated by either theory. Normally, these interme- 
diate flows are encoxmtered when the hypersonic similarity parameter based 
on nose fineness ratio is in the neighborhood of 1. Since this is a range 
of practical interest, additional theoretical methods are needed. 

Some of this need has been fulfilled recently by the hypersonic small- 
disturbance theory (refs. '7 and 8). In its present state of development, 
however, this theory has application only to limited classes of noninclined 
bodies of revolution. Foi- example, due to the series form used to repre- 
sent the pressure distribution, it cannot be applied to the nose-cylinder 
combinations commonly employed for missile bodies. In large part, then, 
the need for a theory applicable at values of the hypersonic similarity 
parameter near 1 still remains. 

The present report develops a theory intended to fiilfill this need. 
This theory is a second-order shock-expansion method. It is developed 
by an i teration procedure which employs the gen erali zed met hod of refer- 
ences JL 'throug h 4 as t ^ first approximation. Expressions are derived 
which define the pressure^qn-nonln’clined— bodies of revolution. Expres- 
sions are also obtained for the normal-force and pitching -inoment deriva- 
tives at zero angle of attack. Predictions of the method, ar.e compared 
with those of other theories and with experimental results. 


SYMBOLS 


A body cross-sectional area 

% 


body base area 
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B 

CijCa 

Cn. 

Cp 

d 

E 

f 


H 

I 

M 

P 

<1 

s,n 


x,r,cp 


X 

a 

3 

7 

n 

A 


fxmction defined by equation (6) 
characteristic coordinates 


normal-force coefficient. 


normal force 


pitching -moment coefficient. 


^o^B 

moment about body vertex 
Q.cA'E^ 


pressure coefficient. 


P - Pq 
^ o 


body diameter 
entropy 
fineness ratio 

(Fineness ratio of the nose section is fn- ) 
total pressure 
body length 
Mach number 
static pressure 
dynamic pressure 

rectangular coordinates (streamline direction and normal to 
streamline direction, respectively) 

cylindrical coordinates (x measured from vertex of body and 
cp from windward meridian) 

center-of -pressure position (measured from body vertex) 
angle of attack 

function defined by equation (12) 
ratio of specific heats (' l.^OO for air) 
flow deflection angle 
fxmction defined by equation (9) 
loading (defined by eq. (l4)) 
function defined by equation (5) 


X 



k 
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Mach angle (arc sine l/M) 

V Prandtl-Meyer expansion angle 

a shock- wave angle 

function defined by equation (13) 

ratio of cross-sectional area of streamtube to that at M = 1 
(see eq. (7)) 

Subscripts 

o free- stream conditions 

1,2, 3, 4 conditions evaluated at various points in flow field 

a afterbody 

c quantities evaluated for cone tangent to the body 

s quantities evalioated by generalized shock-expansion method 

method 

V quantities evaluated at vertex of body 

h quantities evaluated along downstream face of shock wave 

tcv quantities evaluated for cones tangent to body vertex 

tcx quantities evaluated for cones tangent to body at station x 


THEORY 


In the present development of a second-order shock-expansion method, 
attention will be restricted to bodies of revolution. It is recognized, 
however, that the procedure used herein may, in the fut\n*e, find applica- 
tion to other three-dimensional shapes. 

The present method is a refinement of the generalized shock- expans ion 
method of references 1 through k. On the surface of a body of revolution, 
immediately behind a corner, the generalized method represents a first- 
order solution for the flow and the present method gives the second-order 
solution (see Appendix A). Before proceeding, therefore, it is well to 
orient the present analysis with a review of the approximations involved 
in the treatment of the flow about bodies of revolution with the gener- 
alized method. These approximations may be listed as follows (see, e.g, 
ref. 4 ): (l) Distiirbances incident on an oblique shock wave are largely 
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absorbed therein, and hence, reflected distvirbances are negligible; (2) 
the flow appears locally two-dimensional; (3) surface streamlines may be 
taken as meridian lines . In the intermediate range of supersonic speeds 
of interest here, the first approximation is particularly well justified 
(see, ref. S) , and it will not be considered further. As a consequence 
of the second approximation, a solution given by the generalized method 
satisfies the continuity equation only approximately. ^ Although the con- 
tinuity equation does not appear explicitly in the following analysis, it 
is this approximation that is refined by the present meth 9 d. The third 
approximation is one for bodies of revolution only when they are inclined. 
In the present investigation, only bodies near zero lift will be consid- 
ered. Itader this restriction to infinitesimal angles of attack, ein anal- 
ysis has shown that the deviation of true streamlines from meridian lines 
has negligible effect on surface pressures. In the following development, 
therefore, the use of meridian lines as streamlines will be retained. 


Nonlifting Bodies 


The generalized shock-expansion method was developed for nonlifting 
bodies of "revolution from the method of characteristics (ref. 2). This 
development may be summarized with the aid of the equation for the stream- 
wise pressure gradient.^ 


^ ^ ^ 1 Sp 

Ss sin 2[i 3s cos p 3Ci 


In the generali ^j. m ethod the pressure is considered constant along first- 
family Macii’ lines (refs. 1 and !<■) . As a consequence, the right-hand 
member of equation (l) is approximated by zero, and the equation can be 
integrated to yield the well-known Prandtl-Meyer relation. The objective 
of the present analysis is to refine this approximation to the right-hand 
member of equation (l). To this end, consider the flow about a body of 
revolution which has a pointed nose and over which the flow is everywhere 
supersonic. The problem will be simplified by approximating the profile 


^In the treatment of two-dimensional flows, the first approximation is 
used but continuity is exactly satisfied. 

^This equation can be derived directly from the continuity, momentum, 
energy, and state equations with the aid of characteristics theory (see, 
e.g., refs. 2 and 9)* In this form, the equation applies equally well 
for rotational and irrotational flows, requiring only that dE/ds not 
dE/dn be zero. . 
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Of the body with a series of tangents to the original contoxir (see 
sketch (a)). It might be noted that Ferrari (ref. 10 ) suggested a similar 
scheme with a body whose profile was made up of chord lines joining points 



on the original contour. While either approximation is permissible, the 
tangent body was selected here so that the conical flow at the vertex 
will be correct regardless of the degree of approximation used downstream 
of the vertex. 



Sketch (b) 
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For simplification, the derivative, ^p/3Ci . will be approximated with a 
difference equation; thus, along the straight-line element, equation (l) 
may be written (since db/Ss = 0) 

(p) 

. ACiCOs p 

where Ap is the net change in pressure along Mach lines emanating’ from 
the s\irface and ACi is the corresponding length. This equation will be 
solved by an iteration procedinre based on the solution given by the gen- 
eralized method. As previously noted, with the generalized me thod the 
flow is considered two-dimensional and, consequently, no pressu re .change 
is p redi cted along streamlines between the expansion fans .at either end 
of the straight section. While this approximation may be appreciably in 
error for the svirface streamline, it is apparent that the real flow will 
appear more nearly two-dimensional at large distances from the body axis. 

It is reasoned, therefore, that a streamline, well removed from the axis 
(line AB in sketch (b)), can be fo\md along which the pressure will also 
be constant to the accviracy required here..^ For all Mach lines (such as 
CD) emanating from the straight s\n*face then, the pressures at the points 
of intersection with this streamline will be equal. The term, Ap, in 
equation (2) therefore can be written as ki - p, where ki is a constant 
and, of course, p is the varvi ng-.aurface_p r.essirre . The generalized 
method also prescribes that the length (from the surface to streamline ADB) 
and inclinations of all Mach lines will not change when the s\nrface is 
straight. The term, ACiCOs n, therefore can" be represented by a second 
constant, l/kg. Equation (2) thus may be written 


|e = fe(ki - p) 

which can be integrated to yield 

c,p = ki + kae'^f^ (3) 

- V ' 

where ka is the constant of integration. This analysis serves to 
establish the form of the equation representing the pressure distribution 
on an element of the tangent body.'^ It remains now to evaluate the. three 
unknown constants in equation ,(3)» Three known conditions can be employed 
^Examination of characteristic solutions for the flow about cone- 
cylinders indicates that the pressinre variation along streamlines, a moder- 
ate distance from the surface, is markedly less, than that along the surface. 

■*It might be noted that Ehret, Rossow, and Stevens (ref. 11) found 
that pressures on ogives correlated according to the hypersonic similarity 
law could be represented approximately by an exponential function of 
distance. 
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for this purpose. First the pressure, just downstream of the corner, pg, 
can be calculated exactly from the Prandtl-Meyer equations if the pres- 
svire, pi, and the Mach number. Mi, upstream of the comer are known. 
Second, the pressure gradient just downstream of the comer may be cal- 
cialated from the results given in Appendix B. The expression defining 
this 'gradient is 



and Q is the one -dimensional area ratio or 




^ - M 




1 + 


^'M2 


7 + 1 


7+1 

2(7-1) 


( 7 ) 


For the third condition, it is apparent that the pressure on the element 
shown in sketch (b) would approach some limiting value if, rather than 
ending at point 3, the element were extended as indicated by the dashed 
line. If the element were considered to be infinitely long, so as to 
form an extended conical surface, then the only effect the region ahead 
of point 2 will have on the flow at infinity is to form an infinitesimally 
thin layer near the surface across which the entropy varies. It can be 
demonstrated, however, that there is no pressure change through this layer 
and that the flow outside the layer is conical. Consequently, the limiting 
pressure is simply, p^., the pressure on a cone tangent to the original 
body at the same point as the element shown in sketch (b)- (and, of course, 
traveling at the same free -stream Mach number). With these three condi- 
tions, the three unknowns in equation (3) may be evaluated and there is 
obtained . " — 

P - Pc - (Pc - P2)e"’^ 


( 8 ) 
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where { 


^s,4 (Pc - ^2 


(9) 


It is apparent that, in order to apply equation (8), the pressure 
(and Mach n\jmber) on the surface of n oninclined - conss mu sjb_h e kno wn. These 
quantities may be determined from the results of reference 2 or reference 
12. For convenience, the curves shown in figure 1 have been plotted from 
the resvilts of reference 12. 


By application of equation ( 8 ) on successive elements, the pressure 
distribution on the tangent body can be determined. In particular, the 
pressure at each of the points of tangency may be calciilated and appriFd 
to the original body.. The procedure is as follows: First, the elements 

of the tangent body are selected and the coordinates (x.r) of each corner 
determined « The first element is tangent to the body at the vertex, and 
the flow over this element is thus conical. For the first corner, then, 
the pressure, pi, and the Mach nimiber. Mi (see, sketch (b)), are the same 
as^;^the. vertex of the original body. Th e pre ssure, ipa,^ aM the Mach 
nimiber, Ms^, may then be determined with the_ Pj^dtl -Meyer _^imt ions . The 
pressure gradient, (dp/ds)o, may be determined from equation (4) since, 
for the first corner, ( 8 p/ds)i Is zero: The tangent-cone pressure, Pc, 

may be obtained from reference 12 or figinre 1 . With the various factors 
in equation ( 8 ) thus evaluated, the pressure at the tangent point and at 
point 3 (see sketch (b)) can be calculated. In like manner, the pressure 
gradient at point 3 can be determined by differentiation of eqimtion ( 8 ), ‘ 
or 


[tl = (^) (I), <“> 

With the pressure and pres siure gradient at point 3 known (the Mach number 
may also be calculated from the pressure in the usual manner), the factors 
in equation (8) may be determined for the next element. This process is, 
of course, repeated for each element of the tangent body. 


The procedirre just described is not difficult to apply; however, fior- 
ther simplification can be obtained by the use of a tangent 

b ody. This body consists of a cone tangent to the original body at the 
vertex and a conical surface tangent to the body at the station where the 
pressure is to be calculated. With this two-step body, the second S'urface 
is a variable depending on the station in question on the original body. 
For this approximation, equation (8) becomes 




P = Pc - (Pc - Ps)® 


( 11 ) 
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where 


3 = 


X sin - r cos 8y 
r cos ”6 - X sin 6 


( 12 ) 


Bs _ sin 5 a 

(Pc - Ps) sin '5v^j 


(13) 


The subscript, s, denotes quantities at the station on the body as eval- 
uated by the genera lized shnck-fixpansi on-m ethod . With equation (ll) it 
is possible to obtain, very rapidly, a first approximation to the pressure 
distribution. 

The second-order shock-expansion method has been developed to predict 
the pressures on a noninclined body of revolution. In the following sec- 
tion this method will be extended to lifting bodies. 


Lifting Bodies 


For inclined bodies of revolution, a second-order shock-expansion 
method would involve not only a revised expression for the pressures, but, 
in addition, a revised approximation to the shape of the sxarface stream- 
lines. It is recalled from the results of Eggers (ref. 1) that, according 
to the generalized method, surface streamlines may be approximated by 
goedesics. For bodies of revolution. Savin (ref. 3) noted that the per- 
tinent geodesics are simply meridian lines. While this result is exact 
for noninclined bodies of revolution, it is only an approximation in the 
case of inclined bodies. A refined approximation corresponding to a 
second-order method \mdoubtedly could be obtained by graphical integration 
of the momentimi equations employing the pressure distribution given by the 
generalized method. However, it seems at present that this procedixre would 
involve extensive calculations. If attention is restricted to bodies near 
a = 0, it can be demonstrated that the deviation of the true streamlines 
from the meridian lines will not influence siorface pressiores. The approx- 
imation of meridian lines as streamlines can, in effect, be retained and 
relatively simple expressions can be obtained, therefore, for the initial 
slopes of the normal -force and pitching -moment curves. To this end, the 
expression for the normal -force derivative can be written^ 


2jt 

= — / Ar < 


da 


dx 


( 14 ) 


"The subscript, a = 0, has been omitted for simplicity of notation. 
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where A is the nohdimensioneil loading on a thin disk normal to the body 
axis and having vmit radius. This loading A is given by the equation 



<i(p/Po) 

da 


cos cp dcp 


( 15 ) 


The problem then is to evaluate d(p/pQ)/da. The development given pre- 
viously which led to equation ( 3 ) also applies to bodies at infinitesimal 
angles of attack. Equation (8) also applies; however, the variables in 
t his eq uation must be^ oonsidered”as dependent on angle of attack. By 
differentiation of equation (8), there is obtained 


d(p/Po) ,, d(pc/Po) . ^-n d(p2/Po) 

d^ = U - e . ; — + e — 


+ (Pc - Po)e 


-n 

da 


d(p/Pn) d(p 2 /po) 

This equation must satisfy the condition — ^ at t] = 0 

(i.e., X = X 2 ). By the application of this condition to equation (I 6 ), 
the last term (involving dq/da) is eliminated.® The term, d(p 2 /po)/da 
may be evaluated with the aid of the Prandtl-Meyer equation 


d(pi/Po) 

da 


d(Hi) 


Ferri (ref. I 3 ) has shown that the entropy (and hence the total pressure, 

H) on the surface of an inclined cone is constant (independent of cp) . 

When equations (15)» (I 6 ), and (l?) are combined, then, the integrals of. 
the terms involving ■ dHi/da and dH 2 /da will be zero (since J^cos cp dcp=0). 
Equation (15) may therefore be written ® 


(1 - e-"1) 


■n ^ d(p.i/pQ) 

Ai da 


The only terms in equation (I 8 ) that are functions of cp are d(pg/pQ)/da 
and d(pi/pQ)/da. These two terms may be evaluated in terms of the normal- 

dCu 

force derivative of the tangent cone, tcx’ terms of Ai. After 

performance of the necessary manipulations, there is obtained 

©This result indicates that the lifting pressures at small angles of 
attack vary in a manner analogous to that of the pressures at a = 0. 
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A = (1 - e“T)(tan 6) 


dCH 


da 


+ ^ e 


tcx 


(19) 


It is apparent from equation (19) that dCT\r/da for cones must be 
known before the loading A can be evaluated, 
cones are available from 


Fortunately, resvilts for 

have been plotted for convenience 


in figvire 2._ The loading. A, may thus be calculated in the same manner 
as the zero-lift pressures. In this case, Ai, for the first corner is 
dCu' 


simply (tan 6v) 


tcv 


As before, a first approximation to 
two-step body. This approximation gives 


A can be obtained with the 


A = (1- e'l^'l')(tan 6) ^ 

da 


tcx 


+ h. 6. 


V da 


tcv 


( 20 ) 


In .^p endix C , it is shown that eq uation ,(?0) leads to very simple results 
f or~~ ceF£ain common body _ sha pes..=> 

With the loading. A, known, the normal-force derivative may be eval- 
uated by integration of equation (l4) . In like 'manner, the pitching- 
moment derivative can be determined from the equation‘s 


dCm -2n r ^ 

“t — = / Arx dx 

da Agd 


( 21 ) 


A second-order shock-expansion method for bodies of revolution has 
been developed to predict the pressure distribution and the normal-force 
and pitching-moment derivatives at a = 0. The results are relatively 
simple in form and may be applied to a given body with only a moderate 
amount of computations required. Simplified expressions based on an addi- 
tional approximation have also been presented which further reduce the 
amoimt of work required. It should be noted, however, that open-nosed 
bodies and pointed bodies which produce shock waves other than the one at 
the vertex require special forms of the method.® The necessary equations 


■^The contribution to the pitching moment of the variation in local 
axial forces with angle of attack is small for slender bodies (see 
ref. 15) and will be neglected throughout the present analysis. 

®It may also be noted that boattailed bodies present a special problem 

since neither p<, nor -5^ is defined in this cas,e. In practice, how- 

ever, it has been found by comparison with results given in reference I6 

dCNl 

= 2 gives reasonable results for 

tcx 


that the use of 


Pc = Po -5^ __ 

bodies having moderate amounts of boattail . 
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for these cases are contained in Appendix B. In addition, there are 
severeil restrictions on the present method which should be mentioned. 
First, it is apparent that if the exponential variation of the pressures 
is to be valid, then the pressure gradient just downstream of the corner 
must have the same sign as the pressure difference, pc - pa- This condi- 
tion is given by t) > 0 in the ^neral case and by i|r > OTfbf the simpli- 
fied method. There is an additional restriction on the simplified method, 
and that is that the two-step bodies must be real bodies, (i.e., the 
intersection of the two tangent lines must not occvir at negative values 
of X or r). This condition is given by 3 > 0. When n = 0 or 3\|r = 0, 
all equations reduce to those given by the generalized shock-exp^sion 
feunoa. rr.-^ . . 


It remains, of course, to determine the accxnracy of the second-order 
shock-expansion method and to define its range of applicability. There 
are siifficient data available, both from experiment and from character- 
istic solutions, with which the predictions of the method for zero-lift 
pressure distributions can be compared. However, for the case of lifting 
bodies, . sufficient data are not available, and for this reason, the exper- 
iments next discussed were conducted. 


EXPERIMENT 


An experimental program was conducted to determine the initial slopes 
of the normal-force curves and the centers of pressure for a series of 
nose-cylinder combinations. The tests were designed, of course, to permit 
a check on the acciaracy of the predictions of the second-order shock- 
expansion method just developed. It is recalled that the method is 
intended for application at values of the hypersonic similarity parameter, 
MoAn in the neighborhood of 1. The tests cover a range of Mo/fn from 

0.43 to 2.09. 


Appairatus and Tests 


The tests were conducted in the Ames 10- by l4-inch supersonic wind 
txmnel at Mach numbers of 3*00, 4.24, 5*05 and 6.28. For a detailed 
description of this wind tionnel and its aerodynamic characteristics, see 
reference 17. Normal forces and pitching moments for the test models 
were measured with a strain-gage balance. The balance consisted of a model 
support sting on which the moments were meas\n*ed at four points. From 
these four measurements, the normal forces and centers of pressure were 
determined and checked. Measurements were made at nine angles of attack 
from -2° to +4° at each test Mach number. At each angle of attack, the 
values of x/d and Cu/a were calculated. These values were plotted as 
a function of. angle of attack^ and the intercepts at a = 0 of the result- 
ing curves gave the values of dCjj/da and x/d at a = 0. 
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Wind-timnel calibration data (see, ref. 17) were employed in com- 
bination with stagnation-press\are measurements to obtain the stream 
dynamic pressures . Reynolds numbers based on the maximum diameter of the 
models were 


Mach mamber 


Reynolds ntunber, 
million 


3.00 0.79 

4.24 .72 

5.05 .35 

6.28 .15 


Models 


Cones and circiolar-arc tangent ogives of fineness ratios 7, 5, and 
3 were tested alone and with cylindrical afterbodies having lengths of 
2, 4, 6, and 10 diameters. The models were made of polished steel and 
each had a base diameter of 1 inch. 


Accuracy of Test Results 


Stream Mach mmibers in the region of the test bodies did not vary 
more than ±0.03 from the mean values at Mach nxmibers up to 5 •05* A maxi- 
rnimi variation of ±0.05 existed at the highest test Mach number of 6.28. 

The accuracy of the test results is influenced by uncertainties in 
the measxirement of moments and in the determination of the stream dynamic 
pressure and angle of attack. These uncertainties resulted in estimated 
maximum errors in the normal-force derivatives and centers of pressure as 
shown in the following table: . 


Mo 

dCu/da 

x/d 

3.00 

4.24 

5.05 

6.28 

±0.15 

±.15 

±.20 

±.25 

±0.10 

±.10 

±.15 

±.20 


It shovild be noted that, for the most part, the experimental results 
presented herein are in error by less than these estimates. 
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RESULTS AM) DISCUSSION 
Nonlifting Bodies 


The second-order shock-expansion method has developed primarily 

to treat flows characterized by _ val ues. of Mp /fn near unity. Accordingly 
the method has been employed to obtain the. zerb-Tift~pres sure distribu- 
tions at Mp/^n = 1 ^ several different body shapes . ® The resialts are 
shown in figure 3 along with distributions obtained with the generalized 
shock-expansion method (ref. 2). Distributions obtained with the method 
of characteristics (refs. 11, I 8 , and 19 )^ which are considered to be 
exact, are also shown.' It is apparent in figure 3 that the present method 
provides an improvement over the genera,lized method. The differences in 
the distributions obtained with the present method and those obtained 
with the method of characteristics are almost indiscernible. 

In figiore 3(c), comparison is also made with the predictions of the 
hypersonic small-distiarbance theory (ref. 8). The cxnrve shown was cal- 
culated by three terms of a power series representation of the pressinre 
distribution. As noted in reference 8, additional terms will be required 
before this method will accurately predict the pressures on an ogive. 

Even when the additional terms are obtained, however, it seems unlikely 
that the small-disturbance theory will provide a more acciirate estimate 
of the pressures than provided by the present method for the case shown. 
The small-disturbance theory does have a certain advantage in simplicity 
for, if the coefficients of the "series expansion are known, the pressure 
distribution can be calc\lLated very easily. This advantage is partially 
offset by the restriction that the series method requires the body profile 
to have continuous derivatives up to the same order as the nvunber of terms 
used in the series. With this restriction, the theory . cannot be applied 
beyond the nose-cylinder j\mct\ire of the body (fig. 3(c)). 

To investigate the accioracy of the present method at values of Mo/fn 
other than 1, the comparisons shown in figiire 4 have been made. Here, 
the predictions of the present method and those of the generalized method 
axe compared with experimental results for fineness ratio 3 5 tangent 

ogives at Mach niimbers of 3*00> 4.24, and 5-05*^° The values of Mo/fn 
range from O. 6 O to 1.68. The experimental results were taken from refer- 
ence 3* For all cases shown, the predictions of the present method are 
within the accuracy of the experimental data. It is also apparent that 
~ 9In all applications of the present method to curved bodies, the 
tangent bodies employed were formed by elements tan gent to-t he—orlginal 
bodies a!i r~sta ^6hs x/Zn = 0, 0. 1. 0.‘2.'^.'"t T. 1.0. The tangent-body 
approximation is required only if the body profile is curved since for 
cone- cylinders, and for the cylindrical section of any nose-cylinder com- 
bination, the present method yields results in closed form. 

lOpor some of the cases shown in figure 4, the semi empirical methods 
of reference 20 may be used. 
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the predictions of the present method tend to approach those of the gen- 
eralized method as Mg/fn becomes appreciably greater tha n 1 . At 
Mo/fn = 1*68 (fig. ^(f)), for example, the predictions of the two methods 
differ only slightly. 

In figure k, comparison is also made with the second-order potential 
theory (ref. 5) for conditions where this theory is applicable (i.e., 

Mo/fn = 0.60 and O.85). It is somewhat surprising that the present method 
is as acc\jrate as the second-order potential theory even at the relatively 
low value of Mo/fn of O.60. 

The results presented in figures 3 and U indicate that the present 
method fulfills its intended purpose by providing an estimate of the 
pressures on noninclined bodies of revolution for values of Mo/fn near 1. 
At values of Mo/fn as low as O.60 the present method provides res\ilts 
comparable in accuracy with those obtained with the second-order potential 
theory. At values of Mo/fn approaching 2 , the predictions of the present 
method and those of the generalized shock-expansion method differ only 
slightly. It remains now to investigate the applications of the method 
to inclined bodies . 


Lifting Bodies 


The experimental results obtained in the present tests are given in 
tables I and II. Predictions of various theories are also tabulated. 

These include the predictions of the present method (with various approxi- 
mations), the generalized shock-expansion method (ref. 3)^ first-order 
potential theory (refs. 6 and 21 ), Van Dyke's hybrid potential theory 
(ref. 6), and Newtonian impact theory (see, e.g., ref. 22 ).^^ With the 
exception of the two potential theories, all theories have been applied 
throughout the entire range of test variables . The potential theories 
cannot be employed, of course, if the free- stream Mach angle is less than 
the body semi vertex angle. 

Normal-force derivative .- The experimentally determined normal-force 
derivatives and the predictions of the various theories for the bodies 
tested are shown in figirres 5(a) through 5(f). In general, the present 
method predicts the normal -force derivatives at zero angle of attack 

^^Solutions with the second-order potential theory employed in the 
application of the hybrid theory were obtained with the aid of refer- 
ence 23. (Additional results obtained with the first-order and hybrid 
potential theories and with Newtonian impact theory may be found in ref- 
erence 2 k . ) 

^^Curves for the first-order potential theory are not shown in fig- 
ure 5 since, in all except a few cases, the predictions of this theory 
did not differ significantly from those of the hybrid potential theory 
(see tables I and II ). 
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essentially within the accuracy of the data (within about ±0.2) through- 
out the entire range of test variables. In addition, the present method 
appears to provide the most consistently accvirate res-ults of all the 
theories presented in figure 5* The acciaracy of the method at low values 
of Mo/fn can be explained partially by examination of the predictions 
of the method for the limiting case of very slender bodies . In this 
limit, it can be shown from equations (4) and (9) that the term, tj, 
approaches infinity. From equation (l9)^ then, the loading^ A, may be 
written 


A = 2 tan 6 



( 22 ) 


since = 2 (see, fig. 2). With the substitution of this equation 

da tcx 

in equation (l4), there is obtained 


da 




(23) 


This resiilt is, of coiurse, the well-known prediction of slender-body 
theory, which is known to be accurate for slender bodies at low supersonic 
speeds . Thus, the ac c\xrac v of the pres .ent-method-at_low val ues of Mo/fn 
can be attributed, in part, to the fact that it reduces to slen~de ~r=bodv 
theo ry/ ~ln~tEfe~~“rimit I ' ^ 

From the results given in figure several observations can be made 
concerning the accuracy of other theories. For example, it might be 
expected that the potential t heories w ould be more accinrate than the other 
theories when the parameter s/Mq^ - 1 tan by is appreciably less than 1. 
For the fn = 7 cone at Mo = 3 (fig- 5(a-))> however, this parameter is 
only 0.20, and yet, for the longer afterbodies, the hybrid potential theory 
is appreciably more in error than the present method. As found in ref- 
erences 2 through k, the generalized shock-expansion method gives accurate 
results when Mo/f is greater than about 1. Caution should be expressed 
here, however, for the significant parameter is truly Mo/f and not 
Mjj/fn- The results shown in figure 5 indicate that although Mo/fn “^y 
be appreciably greater than 1, for cases where the afterbody is siaffi- 
ciently long to reduce Mo/f below 1, the predictions of the generalized 
method may depart appreciably from the experimental resilLts . In general, 
impact theory gives acceptable results only for nose sections without 
afterbodies . 


Center of pressure . - The experimentally determined centers of pres- 
svire and predictions of the various theories for the bodies tested are 
presented in figure 6. The present method predicts the location of the 
centers of pressure essentially within the accuracy of the data (within 
about ±0.2 body diameters) throughout the entire range of test variables. 
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In addition, the present method again provides the most consistently 
accurate results of all the theories presented. In general, all observa- 
tions made previously regarding the reliability with which the various 
theories predict the normal-force derivatives can also be made in the 
case of the centers of pressure. 

Ranges of applicability .- Several parameters are useful for defining 
the ranges of applicability of the various theories . The ranges of these 
parameters covered by the present tests are shown in the following table: 


Parameter 

Range 

Mo 

3.00 to 6.28 

f 

3 to 17 

fn 

3 to 7 

^a 

0 to 10 

Mo/f 

0.18 to 2.09 

Mo/fn 

0.43 to 2.09 

^/Mo^ - 1 tan by 

0.20 to 2.12 


The second-order shock- expans ion method was foimd to be applicable 
throughout the ranges of variables shown, in the table. Both dCu/da and 
x/d were predicted within ±0.2. The pre sent tests did not reveal the 
limi ts of applicability of the method..— It was indicated , "ho wever , that 
the method may apply to relatively low values of Mo/fn n/Mq^-I tan by) ^ 
since, in the limit of very slender bodies, the method reduces to the 
well-known slender-body theory. Ihe upper limit of the method is dictated 
by the condition specified in the development - namely, T] > 0 (see 

(8))* Calculations have revealed that this condition will be violated 

if tan &y is appreciably greater than 2.5. 

The present tests also reaffirmed the conclusion given in references 1 
through k, that the generalized shock-expansion method is applicable when 
Mq/I is greater than about 1. At values of appreciably greater 

than l,.no significant differences between the predictions of the general- 
ized and second-order methods were found. The ranges of applicability of 
these two methods overlap and thus include most flows about pointed bodies 
of revolution throughout the intermediate- and high-supersonic speed 
ranges. 

Applicatio n of th e potential theories is, of course, limited by the 
condition that JM q ^-1 tan by must be less than 1. Even at the lowest 
values of n/Mq^.-I tan 5y covered by the present tests, however, neither 
the first-order nor the hybrid potential theory was foimd to provide con- 
sistently accinrate predictions of dC^/da or x/d. The calcvilations per- 
formed also revealed no sig nifican t differences in the predictions of the 
two theories at values of n/Mq^-I tan 5y less than about 0.7* 
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Approximations of the Present Method 


As noted in the development of the present method, a simplified solu- 
tion for bodies with curved profiles can be obtained by the use of a two- 
step tangent body. This approximation has been applied to the ogive- 
cylinders of the present tests. By the use of additional approximations 
to the loading. A, the simplified solutions for dCj^g /da _and x/d can be 
obtained, in closed form as discussed in Appendix C. Examples ^ of the 
accuracy of the approximate solutions are shown in figure 7- While the 
approximate methods do not yield results so consistently accurate as those 
obtained with a more complete solution, the approximate methods may still 
be useful to obtain rapid estimates of dCj^j/da and x/d. In this connec- 
tion, these qu^tities can be estimated for ogive- cylinders in a very 
few minutes with the aid of the results given in Appendix C. 

CONCLUSIONS 

A second-order shock- expans ion method applicable to bodies of revolu- 
tion near zero lift has been developed. For noninclined bodies, the pres- 
s\ire distributions obtained with the method were compared with existing 
experimental resiilts and with the predictions of other theories . For 
inclined bodies, the normal-force derivatives and centers of pressure at 
zero angle of attack determined with the method were compared with the 
•predictions of other methods and with experimental results , Cone- and 
ogive-cylinders with fineness ratios from 3 to 17 were tested at Mach 
numbers from 3-00 to 6.28, corresponding to a range of values of the 
hypersonic similarity parameter based on nose fineness ratio (i.e., the 
ratio of free-stream Mach number to nose fineness ratio) from 0.43 to 
2.09* These comparisons led to the following conclusions; 

1. For noninclined bodies, the present method predicts the pressure 
distributions within the accuracy of experimental resilLts . At values of 
the hypersonic similarity parameter based on nose fineness ratio as low 
as 0.6, the present method is as accurate as the second-order potential 
theory. At values of the parameter approaching 2, the predictions of 
the present method differ only slightly from those of the generalized 
shock- expans ion method. 

2. For inclined bodies, the normal-force derivatives and the loca- 
tions of the center of pressure at zero angle of attack predicted with 
the present method are in good agreement with the experimental results 
throughout the entire range of test variables. Within this range, the 
present method yields resiilts more consistently accurate than those of 
other available theories. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 12, I 955 



20 


NACA TN 3527 


APPENDIX A 

POWER SERIES REPRESENTATION OF FLOW ABOUT BODY OF REVOLUTION 


The acciiracy of the present method has been demonstrated by compar- 
isons made over a wide range of flow parameters. It is also informative, 
however, to examine briefly the mathematical accuracy of the method. For 
this purpose, the model shown in sketch (c) is useful. From the vertex 



to point 1, the body is conical. Between points 1 and 2, the surface is 
deflected by a small angle, €. At any point downstream of point 1, the 
physical deviation of the body from a conical surface may be given in 
terras of the angle, e, and the distance, Z\s, measiared from point 1. Simi- 
larly, flow parameters at any point downstream' of point 1 may be expressed 
in terms of e and As. Before developing such an expression, it should 
be noted that for this model (md withi n the restrictio n that the f low is 
everywher e, sunersonic l-.^ the present method provides airi^ac-ir solution for 
the siarface flow in several limits. For example, the present method is 
exact for all values of As when € = p. For As = 0 and As->co, the 
method is exact for all values of For arbitrary values of As and e, 

— of course, the present method is not" exact. However, the general accuracy 
of the method can be demonstrated by expressing flow parameters in the 
form of a Taylor series in the two independent variables, e and As. The 
dependent variable used to define the flow may be any one of several param 
eters. Press\n*e and velocity are among those most commonly used. In the 
present analysis, however, the Prandtl-Meyer angle, v, is considered the 
dependent variable. It should be recognized that the value of the Prandtl 
Meyer angle at a point will define the Mach number, pressure, velocity, 
and other such parameters. We have then the series 
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Each of the derivatives is evaluated at As = e = 0. When € = 0, it is 
apparent that the flow parameters are constant along the surface and 
independent of s. Therefore, all derivatives with respect to. s alone 
are zero. When As = 0, it is also apparent that (^v/d8)i = -1 and that 
all higher derivatives with respect to & alone are zero. We have then ' 
the problem of ^evaluating the cross derivatives. The second-order cross 
derivative, (3 v/Ss ^5), may be evaluated with the aid of equation (Bl4), 
from which (Sv/Ss)2 may be determined; namely. 



It is also apparent that 

\^s ^5 \Ss 


(A2) 


(A3) 


in the limit as 82 ->8x. Hence, by virtue of equation (A2), 



(A4) 


As noted in Appendix B, equation (B14) is not an exact solution for the 
pressure gradient. It can be demonstrated, however, that equation (A2) 
is accinrate to the first order in e' and hence, equation (a 4) is exact. ^ 
^In the derivation of equation (Bl4) a term> 

r 1 /3p \ ds 

2 cos P -2 J A cos p b 

was neglected (see eq. (Bll)).* In the present analysis, both ^p/^Ci 
and the interval of integration, (ss - S4)/b (see sketch (d)), are of 
order e; hence, the neglected term is of order e^. 
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With the substitution of equation (kh) in equation (Al ) , and with the 
application of the other results previously noted, there is obtained 


V = Vi - 6 - 


cos 61 

2ti - 1 


( 'Z - 1 tan 61 ) -1 


(As)(e) + 


0[(As)^(e), (As)(e)^] 


(A5) 


The generalized shock-expansio n method of references 1 through 4 gives 
the resuTlT that \ v = Vi - e.— i Th e general ized method giyes_the Prandtl- 
Mey^ an gle mathe matically accurate to the first order of the independent 
variables e and As and, therefore, iiiimediately downstream of the comer, 
gives "b, first-order solution for the surface flow. The present method 
adds^t^ coefficient of the term involving (As)(e) in equation (A5)_and, 
hence, gives the Prandtl-Meyer angle mathematically accurate to the second 
order of the independent variables e and As. In general, therefore, 
immediately downstream of the corner the present method gives a second- 
order solution for the surface flow, and therefore, it has been termed 
the second-order shock-expansion method. 

The foregoing analysis considered only expanding flows about the 
corner. If € is positive, then the shock wave emanating from the corner 
must be„c 6 nsidered. The result obtained is essentially the same, however." 
For positive e a term of 0(e®) must be added to equation (A5) to account 
for the difference between the Rankine-H\igoniot equations and the Prandtl- 
Meyer equations. Alternately, the term; -e, in equation (A5) can be 
replaced with the change in Prandtl-Meyer angle between points 1 and 2 
as given by the Rankine-Hugoniot equations. The second-order term in 
either case is identical, however, as equation (a 4) may also be obtained 
by differentiating equation (B21). (it may also be obtained by differenti- 
ation of the exact pressure-gradient equation, eq. (B18) . ) 
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APPENDIX B 

EVALUATION OF PRESSURE GRADIENT DOWNSTREAM OF CORNER OF BODY 

OF REVOLUTION 
Convex Corner 


Along a streamline in axially symmetric flow the following relation 
holds (see eq^. (l)) 


^ - Tv - 1 _ -A / ^6 sin sin S' 


Ss 


3^/ cos n ^c\ cos 


From this equation, we may also write 

3p 


^P /^p , SS'X 

SCi Ss/ 


(Bl) 


(B2) 


and 


^6 

aCi 


sin (I sin 5 _ cos u ^ 

r A VSs as 


(B3) 


Consider now the flow in the region of a convex corner on a body of 
revolution as shown in sketch (d). Between points 4 and we may write. 



Sketch (d) 




2 k 


NACA TN 3527 


from equation (B 1 ) 


'P5 


^ - (s= - 64) - r -A_ 

A J A cos [X VSCi/ 


ds (bU) 


If points 1 + and 5 are near to the surface, equation (b 4 ) may be approxi- 
mated by 


/ ^ 
/ A 


Pc; - P2 P4 - Pi 


- (S5 - 6. 


Pi 


^ J A cos |i \ 3 Ci/ 


ds 


(B 5 ) 


Since the flow between points 1 and 2 is strictly of the Praudtl-Meyer. 
type, 





Pi 


iviBr •ys '-r ^2 - 83 




(b6) 


We may also write from equations (B 2 ) and (B 3 ) 


Ps - Pa ^ 


Aa As V^Ci 


P4 - Pi _ /'^P \ a = 


Ai ^ Ai v^Ciy Ai 


8 . ■ 


S. ■ h 8.. 


(B 7 ) 

(B8) 


65-62 + ( ) b = 62 
yoCiy 


sin nasin 62 cos 112 


.r. 


b - 




^ - A; 






(B 9 ) 


R _ R /86\ 

04 - 61 + ( J a = 


sin p^xsin 6x pos ^x 

5 i -p, a 




Hi 


IJ 


a 


(BIO) 
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When equations (b 6) throiagh (BIO) are substituted into equations (B5)> 
there is obtained 


2 cos qg Y^\ 
A2 





sin liasin 62 
r 


b 


sin ij-xsin 5i 
_ 


a = 



1 

A cos M. 



ds 

(Bll) 


If it is assumed that a first approximation to the flow is given by the 
generalized shock- expans ion method, then the right-hand member of equa- 
tinn t'BTl') ma v— be-negleet-edv Equation (Bll) may thus be written""" 



A2 

2r cos |i 2 



I 

sin |j.isin 6i 



sin [i 2 sin 


+ 


cos \lx 
cos P 2 Yv 



(B 12 ) 


In the limit as the streamline between points i)- and 5 approaches the 
sixrface, the ratio, a/b, may be evaluated in terms of the one-dimensional 
area ratio 


a _ /sin |j.2\ ^ 
b ~ \sin niy ^2 


(BI3) 


With the substitution of equation (B13) into (B12), there is obtained 
after combination of terms 
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'P 


H^j - [jr^ 

\ /2 \ ^ 


* I [d 


jp\ _ -V Y66 

v6s / 6s 


(bi4) 


where 


B = 


7pM‘ 


|2 , 


2(vr - 1 ) 


(B15) 


and, of course 


A.= — 


2rp 


sin 2 |i 


(BI 6 ) 



1 

M 




7 + 1' 


(y+i) 

2(7-1) 


(B17) 


Equation (B14) represents only an approximate evaluation of the 
pressure gradient. More exact evaluations may be found in references 10 
and 25 . These more exact results, of course, require m;mierical or graph- 
ical integration. 


Concave Corner 


In most cases, the tangent bodies used in the application of the 
present method will have convex corners. There is a possibility that 
concave corners may be encountered. In the event that the original body 
does not have sharp concave corners, equation (Bl^) will still suffice 
since the flow along the surface is still isentropic. However, if the 
original body does have sharp concave corners, then the pressure gradient 
for this case will also be required. This result can be obtained in the 
same way as equation (B14); however, the shock wave emanating from the 
corner must be considered after the manner described in reference 9 . The 
expression defining the pressure gradient in this case is 


NACA TN 3527 


27 



[(7 + l)tan( 6 (i - 6^)cos(gt - &q) - sin(a - &u)]Mu sin^(a - + sin(a - 6 ^) 

+ [1 - 2 sin^(a - 5u) + 2 tan(6d - 6u)sln(a - 6u)cos(a - 6u)]Mu^sin^(cr - 6u) 
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In these equations, a is the shock-wave angle with respect to the body 
axis, and (SH/Sn)^ is the variation of the total pressure normal to the 
surface just upstream of the shock wave. The subscript, u, refers to 
conditions upstream of the shock wave, and the subscript, d, refers to 
conditions downstream. Equation (BI 8 ) represents the exact solution in 
the usual sense . All effects of the interaction between shock waves and 
Mach waves are therefore included. In order to be consistent with other 
parts of this analysis, these effects should be neglected. In addition, 
since equation (BI 8 ) is intended for application to a tangent body, the 
body curvatures, ( 8 &/ 8 s),ji and ( 86 / 8 s)<i, will be zero. It may also be 
noted that the first step- of the tangent body is a cone tangent to the 
vertex of the original body. For this approximation then,' there will be 
a small layer near the surface of the tangent body for which ( 8 H/ 8 n) = 0 
With these approximations applied to equation (B18), a simplified result 
can be obtained which will suffice for the present purposes. 



r tan 

_ ^ 

sin(g - 5u)sin 5u 

V^s/cl 

tan(a - 5d). 

r 

sin(a - 5d) , , J 



Bd sin(g - 5u) 
^[Bu sin(a - 6 (i) 


+ 




cos(a - 5u)tan |i(j^ 
sin(a - 6 ( 3 ^) 


(B21) 


For a body with a concave corner, a special form must also be used 
loading. Just downstream of the corner and before the first 
convex^ corner 
/ 


I A = (1 - e"^)tan S 
I 




tcx \Pu 





(B22) 


and, thereafter. 




A = (1 - e"^)tan 5 


dC. 


N 




da 




'^i + (rrr"~Tir) ^^u 

tcx ^ 


t>2 A 2 Pi 

Po Ai Bo 


where 


G = 


(Pd/Po) 


[Mu^sin^(a - &u) - l]‘ 


[(7 - 1)Mu sin^(a - &u) + 2]Mu sin^(g - 5u) 


(B23)’, 


(B2k) 
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The corresponding equations for open-nosed bodies of revolution are 
similar. The pressure gradient at the leading edge may be determined 
from equation (B21) with Mu = Mq, &u = (^/^s)u = 0. The loading 

on the exterior surface is given by 


A = (1 - e”’l)tan S 


da 


A2 

+ T" ® 

tcx 


'T. /P2 ^ 

■ (po " Ai Po 


JA-v 


(B25) 


where Ay is the loading at the leading edge, or 


Av = 


4 sin CyCOS Oy 


(r 


r sin(ay - 6y)cos(0y - &v) ^ 

f ^ 

^ COS^(Oy - 6y) 

L sin Cvcos ay 

^ + V 

Mo^sin^Cy 


(B26) 


and J is defined by 


(Mo^sin^Oy - 1)" 


J = 


(Pv/Po)^o^®i^®'^v[(A - l)Mo^sin2ay + 2] 


(B27) 


For bodies with concave corners, and for open-nosed bodies, the total 
pressiore is not constant on the surface when. the bodies are inclined. 

Th is va riat ion in surface total pressure leads to the term involving Aq 
in^q_uati^n (B23;) -and^'the term inwlving Av in equation (B25). 
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APPENDIX C 

EXTENSIONS OF THE APPROXIMATE METHOD 


This analysis is based on the approximate or tvo- stan p t-P- 

viously mentioned. The basic equations of this method are equations (ll), 
(12), (13)> and (20). Before proceeding with this analysis, it is con- 
venient to write down the expressions for the fvinction p (see eq. (l2)) 
for several types of bodies. These expressions are presented in the 
following table : 


Expressions for p 


Body 

For nose section 

For cylindrical afterbody 

Any body 

X sin 6y - r cos by 

2fnsin ^v(^J -- cos 6y 

r cos 6 - X sin 6 

Cone- cylinder 

(Not required) 

- 0 

nTi + 

Tangent- ogive - 
cylinder 

1 


. 1 + ^fn^ 

Tangent -par ab olo id 
cylinder 

H ^ 


y fn^ + 1 

nT f n^ + 1 


In general, the equations for the normal-force and pitching-moment 
derivatives may be integrated in two parts - one part for the nose sec- 
tion and one for the afterbody. Thus, with the loading defined by 
equation (20) 


da 


da 


+ Gie 


-Ggfn 


nose 


0 



(Cl) 



da J 


dCm 


da 


nose 


^ e ^j^(l + Gsfn) - (1 + Gafn + Gafa)® 


(C2) 
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where 


Gi = 


2 Psa sin 2|iy 
'Ifacos 6y Pv sin 2|isa^ 


do/ -tcv 


(C3) 


and 


G2 = 2 \lfaSin 6 y • (C 4 ) 


The additional subscript, a, refers to fimctions evaluated for the after- 
body (i.e., 6 = O). Thus, from equation (13)^ 


y(Psa/Po)Msa^ fiy 

2(1 - Psa/Po)(Msa^ - 1) “sa 


The terms Gi and Ga are fimctions of Mq and 5y alone. These 
functions have been evaluated and the results are shown in figure 8. For 
the special case of cone-cylinders, equations (Cl) and (C2) represent a 
closed solution of the general method as well. 

By the use of an additional approximation to A, results in closed 
form can also be obtained for ogival nose sections. Such an approximation 
is 


A = Gatan 5- 


dC 


N 


do, 


+ G4tan 5 + 


[d- 


G3)fS 


tcv 


(tan 5) 


tcv 


tan Sv 


(C6) 


where 


^3 


Pv 


sin 2qy 
sin 2nsa 


(C7) 


and 


G4 = 2(1 - e‘^n) 


(C 8 ) 


When equation (C6) is substituted in equations (l4) and (21), equations 
are obtained in closed form for dCw/da and dCm/da. These equations 
involve constants which are complicated functions of the nose angle Sy 
(or nose fineness ratio ^ i*n)* These functions can be expanded in a series 
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in terms of 6 y; the leading terms of these series are constants independ- 
ent of 6 v. In view of the approximate nature of this analysis, the use of 
the leading terms will suffice. Thus there is obtained 


. 

da I ogive 



+ 4 (1 + ^3) 



(C9) 


1 



(3 + 22G3) 


doc 

ogive 

■ 15 " 

15 

V^-Acv 


(CIO) 


To^_; ^e accvuracy of this analysis , these equatio ns also r epresent the solu- 
tions_ f.or--.a... tangent . paraboloid These equations have been "evaluated for 
a range of Mach numbers and nose fineness ratios . The resiilts are pre- 
sented in figvire 9. It is apparent that with the aid of equations (Cl), 
(C 2 ), (C 9 )^ (C 10 )> and figures 8 and 9 , dC^/da and dCm/da for ogive- 
cylinders can be evaluated approximately in a few minutes . 
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F/gure 2.- Normal-force derivative at a for cones (ref. 14). 
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Figure 3- Surface-pressure distributions for various bodies at 
Mo^fn - t ond a - O'*. 
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Afterbody fineness ratio, fg 


(a) Cone, f,,- 7 


Figure 5- Normal -force derivatives for various cone- and ogive-cylinders 
at Mach numbers from 3.00 to 6.26. 
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Figure 5- Continued. 
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Figure 5- Continued. 
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Figure 5.- Continued. 
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Figure 5.- Concluded. 
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Figure 6- Centers of pressure at. a -O** for various cone- and ogive- 
cylinders at Mach numbers from 3.00 to 6.28. 
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Figure 6.- Continued. 
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Figure 6- Continued. 
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Figure 6- Continued. 
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Figure 9.- Curves defining normal -force and pitching -moment derivatives 
for tangent ogives as predicted by approximate dosed solution of 
Appendix C. 
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